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Abstract. Exponentially localized solutions of the Klein-Gordon equation for two 
and three space variables are presented. The solutions depend on four free parameters. 
For some relations between the parameters, the solutions describe wave packets 
filled with oscillations whose amplitudes decrease in the Gaussian way with distance 
from a point running with group velocity along a straight line. The solutions are 
constructed using exact complex solutions of the eikonal equation and may be regarded 
iQ} ' as ray solutions with amplitudes involving one term. It is also shown that the 

i-C . multidimensional nonlinear Klein-Gordon equation can be reduced to an ordinary 

differential equation with respect to the complex eikonal. 

Translated from Zapiski Nauchnykh Seminarov POMI, Vol. 275, 2001, pp. 187-198. 
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Introduction. The construction of various highly localized solutions of the wave 
equation 

<j>U - 4>xx - 4>yy - 4>zz = (1) 

is a subject of several publications (see [1-8]). In particular, in [1] a solution that is 
exponentially localized in the vicinity of a point running with velocity of light is given. 
In this paper, we generalize the result of [6] by presenting a family of localized solutions 
that includes the one given in [6] as a special case. 

Based on the results obtained for the wave equation, we construct a family of 
particle-like solutions of the Klein-Gordon equation 

h 2 (uu — u xx ~ u yy) + u = 0, h = const. (2) 

These solutions have finite energy and describe wave packets with central frequency uo 
and wave number k, where uj 2 = k 2 + 1/h 2 . Their amplitudes decrease exponentially 
with distance from a point running along a straight line with group velocity v = du/dk. 
By analogy with the solutions of the wave equation, we call them Gaussian wave packets. 

All the solutions of the Klein-Gordon equation from this class can be represented 
in the form 

u = Af(iS/h), 

where the function S satisfies the eikonal (or Hamilton- Jacobi) equation 

S? - S 2 X - S 2 = 1, (3) 
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the amplitude factor A does not depend on the coordinates, S and A are independent 
on h, and / is expressed in terms of the Hankel function. 

The comparison of our results with those available in the literature shows that 
one of the solutions that we found for the three-dimensional Klein-Gordon equation 
coincides with the solution given in [7] . We also consider in detail another solution from 
the constructed family, which depends only on one variable S that is one of the exact 
complex solutions of the eikonal equation. It turned out that the search for a solution 
dependent only on S of the nonlinear Klein-Gordon equation is reduced to the solution of 
an ordinary differential equation. This is also true for nonlinear Klein-Gordon equations 
with arbitrary number of space variables. 

In constructing particle-like solutions of the wave equation, we followed the idea 
by Ziolkowski [3] and sought such solutions in the form of a superposition of Gaussian 
beams, which are solutions localized near rays. The latter solutions were found for the 
first time in papers of Brittingham [1] and Kiselev [2] and belong to the class of relatively 
undistorted waves, in the terminology of Courant and Hilbert [10]. The construction of 
the solution of the Klein-Gordon equation employs the simple observation that, taking 
the Fourier transform of a solution of wave equation (1), say with respect to z, 

/oo 
dz(p(x,y,z,t)e iz/h , (4) 
-oo 

we obtain a solution of the Klein-Gordon equation (2). 

Generalized Gaussian packets for the wave equation. We start from the solution 
of wave equation (1), described by the formula 
, / t \ exp{ig0i} 

(p b (x, y, z, t, q) = -rz . . , . ( 5 ) 

{(3 - is-i) 1 ' 2 ^ - ie 2 ) l/2 

(see [1, 2, 6]), where the notation 

2 2 

e 1 = x-t + -^— + -^—, p = x + t, (6) 

p — isi p — %e 2 

is introduced. The function 0& satisfies (1) for any q, e±, and e 2 , and in the case of 
£i > 0, e 2 > 0, and q > it is a Gaussian beam, which means that it is localized in the 
Gaussian way in the vicinity of the x axis. We seek particle-like solutions of Eq. (1) in 
the form of a superposition of Gaussian beams: 

POO 

<j>M{x,y,z,t)= / dqF("\q)Mx,y,z,t,q),(7) (7) 
Jo 

where F^(q) is a particular function depending on the parameter v. We put 

F <yV \q) = aq-"- 1 e- e( - q+a3/q) , (8) (8) 

where u, a, and e are arbitrary constants, a > 0, e > 0, and a = (4ecj 2 ) !y /(2y / 7r). It can 
easily be shown that (7) is reduced to an integral representation of the Hankel function 
of the first kind [12] and 



^^..^(l-^i) , (9) 
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where C = i2 v ~ 1 ^pn. It is worth noting that s satisfies the Hamilton- Jacobi equation 
s t = s l + s l + s 1 f° r w ave equation (1). We note that for v — 1/2, formula (9) yields a 
solution of the wave equation presented earlier in [6]: 

exp \ —2oe\/\ — iQi/e \ 
0(l/2) p (x, y, z, t) = I , > (10) 

- i£i){(3 - is 2 ) 

This solution depends on four free parameters e, £i, e 2 , and c- It is established in [6] 
that if all these parameters are positive, it is localized in the Gaussian way near the 
point x = y = and z = ct that runs with velocity of light c = 1 along the x axis. The 
asymptotics of the solutions of (1) of the form (9) with respect to the large argument 
have the same exponential factor as (10). Arguments similar to those adduced in [6] 
prove their localization. Therefore, (9) is a localized solution generalizing (10). 

Gaussian beams for the Klein-Gordon equation. Here we give a solution u b of 
the Klein-Gordon equation, which has a Gaussian localization near a ray, in order to 
use it in construction of particle-like solutions of the Klein-Gordon equation. To find 
such a u b , we calculate the Fourier transform (4) with respect to z of expression (5): 

x ^e iw ^- £2 ^ 4 ^exp{iQq-i(3/(4qh 2 )} . . 

u b {x, y, t; q) = , , (11) 

y/q VP - IE X 

where 

V 2 

e = x-t+ - — . (12) 

p — %E\ 

This solution was found first in [8]. We call such solutions of the Klein-Gordon equation 
Gaussian beams, by analogy with solutions of wave equations. 

Gaussian packets for the Klein-Gordon equation. Taking the Fourier transform 
(4) with respect to z of the both sides of Eq. (7), we obtain solutions of the Klein-Gordon 
equation in the form of an expansion in the beam solutions u b from (11): 

pOO 

uf\x,y,t)= dqF^(q)u b (x,y,t;q). (13) 
Jo 

Substituting u b of the form (11) and F^ v \q) of the form (8) into (13), we arrive at 

M _ Sp +1/2 Hl% /2 (S p /h) 

< (x,y,t) - d ^ _ — jj^-p _ ^ 2h2 - — y^j- 2 (14) 

with 

S p = i [(& + ie)({3 - 4iea 2 h 2 - ie 2 )] 1 /2. (15) 

The function S p does not depend on v. Formula (14) yields a family of exact solutions 
of the Klein-Gordon equation, which are all particle-like as will be seen in what follows. 
This family depends on four free parameters v, e, e±, and a. In what follows we use a 
different and simpler parametrization of (14). The constant factor in (14) is equal to 
d = n(8ea 2 h 2 + 2e 2 y- 1 / 2 h u+ ^ 2 e l ^ u+i yV2. 
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The ray interpretation of Gaussian beams and Gaussian packets. The 

introduction of a new free parameter k by the relation k = qh enables us to rewrite 
(11) in the form of a ray series, which is reduced to only one term 

exp{iS b /h} 

U b = C bJ75 : M/o » °b = COnst, (16) 

with 

S b = nQ-^-. (17) 

Here, the complex phase function S b is independent of h and satisfies eikonal equation 
(3). We also note that u b may be regarded as a reference solution for the asymptotic 
construction presented by Maslov in [11]. 

The solutions u p ^ can conveniently be written in terms of a new free parameter 
7 = Aea 2 h 2 + e 2 instead of a and e 2 and /i — v + 1/2 instead of v. Now, 

The complex phase function 

S > p = 2[(0 + ^)(/5-2 7 )] 1/2 (19) 

is also independent of h and satisfies eikonal equation (3). In the case of a half-integer 
/i, the Hankel functions are reduced to elementary functions. For the values /i — 1/2 
and \i = —1/2, we have 



«i 0) = c 2 



exp {iSp/ h} 



and 



„(-, = «p{^/ft} /|zji (20) 

respectively, where C 2 and C 3 are constants. These solutions, as well as (16), are ray 
solutions reduced to their zero-order term. 

Properties of the solutions. First we show that the solution u b in (16) is localized 
near the x axis. Separating the real and imaginary parts in the complex phase function 
S b from (17) and introducing the notation 



«=^ + T, A, = ^-£±^,(21) (2!) 

we obtain 

i5 6 = i(«x-a;*)-^ + ^-7-- (22) 

For \(3\ £i, the solution u b describes a wave with frequency uj = u/h and wave 
number k = k/h, which propagates along the x axis and decays in the Gaussian way 
with distance from the axis. The degree of localization near the x axis is determined by 
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the parameter hA y . The solution shows the Gaussian localization near the x axis for 
any finite values of f3, and the degree of localization decreases as j3 grows. 

Now we describe the asymptotic behavior of u p ^ l ^ 2 \ As has already been 
mentioned, all the functions S p are independent of fi. The asymptotic expressions of 
the Hankel function for large values of \S p \/h differ only in the phase factor, and all 
solutions of our family behave as 

„(/*-i/2)/ f \ _ c S p 1/2 exp{tS p /h} 



with C 3 = d^lj^e-^/ 2 -™^. 

Now we investigate the dependence of S p on the space variables and time. 
Separating the real and imaginary parts in the subradical expression from (19), we 
obtain 

1/2 



iS p — 



(24) 



where /c and Cj are defined by (21) with k = ^y/Ae. We now prove a rough estimate 
valid for any value of t, which implies that any solution from the family described 
by (14) has finite energy. First we note the simple fact that for real a\ and a 2 , 

/ \l/2 

ReVai + i&2 = ( (ai + y a? + a|) /2 ) , whence 



Re\/ ai + ia 2 > ai for ai > 



and thus 

?2 



|Re(^ p )| > 



_> ,2 2 /3 2 + 7£i llJ 



7^ + x - r + y 



(25) 



/5 2 + e? 

Estimate (25) shows that for a fixed time and sufficiently large values of the space 
coordinates, the asymptotic expression for the Hankel functions is applicable and the 
behavior of solution (18) is described by formula (23). Therefore, the absolute value 
of any solution |wp' u_1//2 - ) | decreases exponentially with a rise in the coordinates, and 
therefore all these solutions have finite energy. This estimate is valid for all values of 
the space coordinates in case of \t\ < 7s. 

Now we prove that, under certain restrictions on x, y, and t and some conditions on 
their parameters, the solutions Up^ 1 ^ 2 ^ describe wave packets localized in the Gaussian 
way. We expand iS p in a series, assuming that all subradical terms dependent on x, y, 
and t are small in comparison with £7. Keeping only the terms linear and quadratic in 
x, y, and t, we find 



iS p ~ i(kx - ujt) - y/^i — 2 --, (26) 



(x — vt) 2 y 2 
A 2 "A 

where v stands for the group velocity v = du/dk = k/u and 
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Figure 1. Real part of in conventional units: t = (left), t = 10 (right). 



We recall that k = y/^y/Ae. 

Under the condition y^fe = 2ke ^> h, the asymptotic form of the Hankel function 
can be used in (18) for all values of x, y, and t. Formula (26) can be applied for not-too- 
large x, y, and t. However, for some values of the parameters k, e, and E\, the domain 
of its applicability can be larger that the widths VhA x and \fhA y of the packet. In this 
domain, all the solutions from the family tip 1-1 ^ demonstrate Gaussian localization 
near a point running along the x axis with group velocity v. Such a behavior may occur 
for a certain relation between all the parameters, which we write below in an interesting 
special case. 

On the Fig. 1. the real part of u p °^ given by (20) is presented in conventional units: 
t — (left), t — 10 (right); in both cases, h = 0.3, e = 10, and k — 1. 

A simple example of a Gaussian packet. An interesting solution of the Klein- 
Gordon equation (2) is obtained from (1) by setting £1 = 7 = 4/t 2 e. The solution u p ~^ 
depends only on one variable S p , and, up to a constant factor, we have 

J,Sp/h 



U 



= a, 



(28) 



In this case, the phase S p is given by 

S p = i^b 2 + x 2 + y 2 -t 2 + 2ib(ut - kx) 

= i^/(x-ibk) 2 + y 2 - (t - ibCo) 2 (29) 

with b = ^Ffs = 2ks. Here, k can be expressed in terms of k, see (21). However, it is 
more convenient to regard k and b as independent parameters. Then, u> = \Jk 2 + 1. 

The expansion of S p takes the form (26) with A x = y/b/u and A y = \p2h. We 
describe the conditions for its validity. The expansion of (29) is justified for 1 « 6, 
y b, t <^ b, t <ti b/u , and x b/k. Since Q > 1, the conditions t b/u and 
x mm(b,b/k) are more restrictive. It can be seen from (29) that the solution is 
localized in the Gaussian way near a point running with group velocity if the transverse 
and longitudinal widths VhA x and \^hA y , as well as the time width y/hA x /v, are 
all far less that the domain where the root in (29) can be expanded, namely, for 
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\fhhuj <C min(6, b/R), \Jlhh <C b, and \ /r hb/R <C b/uo. Here, the most restrictive 
condition is 

h < b%. (30) 

The case of three space variables. All of the above considerations can be generalized 
to three or more space variables. Here we present the results for the three-dimensional 
case and compare them with the results found in [7]. We turn to the equation 

h 2 (utt — u xx — u yy — u zz ) + u — 0, h — const. (31) 

The construction of its localized solutions is based on the expressions for localized 
solutions of the four dimensional wave equation 

4>tt - 4>xx - 4>yy ~ 4>zz ~ 4>z lZ1 = 0. (32) 

Axially symmetric Gaussian beams for this equation are described by the expression 
, / i \ exp{ig0!} 

fe (x, y, z, zi, t, q) = — —z . 33 

In contrast to (6), here 

2 2 2 

Qi=x-t + -V— + -^— + -t^— , P = x + t. (34) 
p — is\ p — ie 1 p — ie 2 

We construct generalized Gaussian packets for Eq. (32) by a formula similar to (7), 
starting from the Gaussian beams and using weight function (8) in the same way as we 
did this previously. Then we take the Fourier transform 

/oo 
dzx4>{x,y,z,z x ,t)c lZlth (35) 
-oo 

of the both sides of a formula analogous to (7). In doing so, we find solutions of the 
Klein-Gordon equation (31) describing the Gaussian particles 

u (^-i/2) _ c S P H P( S p/ h ) m) 

where C4 = const and S p is defined by (19) in which 

e = x-t+-^- + -^—. (37) 

P — IS\ p — IS\ 

The Gaussian particles u^ 1 ^^ in the 3D case differ from the corresponding solutions 
in the 2D case (18) by a power of (5 — isi in the denominator and by the additional term 
involving the variable z in phase (37). The solution 

Up - ° 5 (J3 - ie 1 )S p P 1 p/ 1 " 5 (J3 - ie^e + ie)W { } 

with C5 = const was found earlier in [7]. This is a Gaussian particle from (36) for 
(j, = —1/2. In the three-dimensional case, as well as in the two-dimensional case, a 
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solution dependent only on the variables S p exists. This is u v 3 ^ 2 \ where we must set 
Si = 7: 

< 3/2) = Ce H ' ){ i Sp/h \ (39) 

Op 

where C 6 = const. 

On nonlinear Klein-Gordon equations. The above construction can be used in 
solving nonlinear Klein- Gordon equations 

h 2 (u a - u xx - u yy ) + <p{u) = (40) 

in the two-dimensional case and 

h 2 (u tt - u xx - u yy - u zz ) + (f>(u) = (41) 

in the three-dimensional case. Seeking solutions in the form u = F(s), where s = S p /h, 
we arrive at nonlinear differential equations 

F ss + -F s + 0(F) = (42) 
s 

and 

F ss + -F a + 0(F) = (43) 
s 

in the case of (40) and (41), respectively. 
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